First-principles study for the adsorption of segments of BPA-PC on a-Al2O3(0001) 
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We have studied the adsorption of bisphenol-A-polycarbonate (BPA-PC) on the a-Al2O3(0001) surface 
using density-functional theory (DFT) with van der Waals (vdW) corrections. The BPA-PC polymer can 
be divided into its chemical fragments which are phenylene, carbonate and isopropylidene groups. We 
have calculated the adsorption energy and geometry of the BPA-PC segments that consist of two to 
three adjacent groups of the polymer. Our DFT results show that the adsorption is dominated by the 
vdW interaction. It is also important to include the interaction of nearest-neighbor groups in order to 
provide a realistic environment for the adsorption of the polymer onto the surface. Our results also show 
that the BPA-PC molecule attaches to the alumina surface via the carbonate group located in the middle 
of the molecule chain. 

PACS numbers: 68.43.Fg, 68.47.Gh 



The bisphenol-A-polycarbonate (BPA-PC) molecule 
studied in this work is an important industrial polymer 
used in different composite materials together with, e.g., 
metals. The molecule has extensively been studied both 
experimentally ~ and computationally ~ . The interac- 
tion of BPA-PC or its chemical groups, phenylene, car- 
bonate, and isopropylidene, with crystal surfaces has also 
been studied, and they provide a good model system for 
studying the properties of surfaces and interfaces 8 ' 11-18 . 

Density-functional theory (DFT) has become an impor- 
tant tool for the study of matter at the electronic-structure 
level, including processes on surfaces such as the adsorp- 
tion of atoms or molecules. Even though the covalent in- 
teraction is crucial for the intramolecular binding or the 
binding of single atoms on surfaces, large molecules with 
closed electron shells are interacting relatively weakly 
with the surfaces. Large, especially organic, molecules 
interact with the surface mostly via the van der Waals 
(vdW) interaction. This electric dipole-dipole interaction 
has a longer range compared to the covalent interaction, 
allowing the molecule to interact with the surface compar- 
atively far from it ' . However, the conventional DFT 
includes only the covalent interactions and thus it fails 
to give the correct adsorption energy of large molecules 
on surfaces. Recent developments have made it possible 
to include the effects of vdW interactions in the DFT de- 
scription of the large molecule systems. 

In this work we extend a subset of our previous DFT 
calculations for the groups of BPA-PC (Ref. ) by taking 
into account the interactions of nearest-neighbor groups 
as well as the vdW interactions. Both of these factors 
can be crucial in understanding the adsorption processes 
of large molecules, with macromolecules (polymers) being 
an extreme example. To show this, we have calculated 
the adsorption energy and geometry of different segments 
consisting of two to three adjacent groups of the BPA-PC 
polymer chain on the a-alumina (0001) surface. 

In the numerical calculations we combine two differ- 
ent DFT approaches depending on the purpose described 
below. Compared to adsorption calculations for small 
molecules or individual atoms, adsorption calculations for 
large molecules offer a number of challenges. First, in or- 



der to fit the molecule into the simulation box such that 
it does not interact with its own mirror image due to peri- 
odic boundary conditions, the box has to be large. Second, 
large molecules, especially segments of polymers stud- 
ied here, have large conformational degrees of freedom, 
meaning that relaxing them requires a lot of time. These 
two issues, combined with the fact that one needs a simi- 
lar level of accuracy as for normal adsorption calculations, 
essentially makes a traditional plane-wave DFT approach 
unfeasible. In order to perform these calculations, we 
have to make quick scans with a relative fast linear com- 
bination of atomic orbitals (LCAO) DFT method first and, 
after that, refining the interesting cases with a more ac- 
curate but significantly more expensive real-space finite- 
difference DFT method. This combining of two codes is 
significantly easier by controlling both of them via a uni- 
fied high-level interface 21 . 

For the LCAO calculations, we used the Siesta code , 
using single-f polarized (SZP) orbitals and a single k- 
point at the T point. The many-body effects were approx- 
imated with the revised version of the Perdew-Burke- 
Ernzerhof ' (RPBE) form of the generalized-gradient 
approximation (GGA). The atomic coordinates were re- 
laxed until the maximum force was less than 0.04 eV/A. 
The finite-difference calculations were done using pro- 
jector augmented waves (PAW) as implemented in the 
real-space finite-difference code GPAW 21 ' 27 ' 28 . We used 
a grid spacing of approximately 0.19 A, and a 2x1x1 k- 
point mesh using the Monkhorst-Pack scheme ' . The 
exchange-correlation and vdW interactions were taken 
into account with a revised version of the PBE GGA ;4 ' 25 ' 31 
(revPBE), and a self-consistent implementation of a fast 
Fourier transform (FFT) based method based on the 
vdW-functional of Dion et al. 33 ' 34 . The systems were re- 
laxed until the maximum force on all atoms were less 
than 0.02 eV/A. The Fermi level was smeared using Fermi 
smearing with a smearing parameter of 0.1 eV. 

The a-Al 2 O 3 (0001) supercell contained 96 Al and 144 O 
atoms in 12 Al and 6 O layers forming an 8x4 hexagonal 
surface supercell with Cartesian dimensions 8.3 Ax 19.1 
A x 27.0 A. The height of the slab was 12.1 A, and there 
was a 14.9 A vacuum region. The lowest Al and O layers 
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Figure 1. Structural formula of the BPA-PC monomer. The 
monomer consists of the groups P, C, P, and I, respectively. 



were fixed, while the rest of the slab was relaxed. The 
surface was Al-terminated, as previous calculations have 
shown this to be the most stable termination , although 
there are some conflicting experimental results suggest- 
ing a preference for a mixture of different terminations . 
For some calculations, depending on the adsorbate config- 
uration on the surface, the cell size was doubled to 16.6 A 
x 19.1 Ax 27.0 A. 

A BPA-PC monomer can be seen in Fig. 1. The 
monomer consists of three different chemical groups, 
phenylene (marked with P, chemical formula C6H4), car- 
bonate (C, CO3), and isopropylidene (I, C3H6). As the 
BPA-PC polymer chains are usually terminated at both 
ends by a P group , for the adsorption a PC segment 
was simulated. The adsorption of the PI segment has not 
been simulated, as the I group sterically prevents the P 
group from getting close to the surface. An IC segment 
has not been studied either, because it does not contain 
any P group. For the intra-chain adsorptions, segments 
consisting of three groups were simulated, but the chains 
of length of four groups or longer have not been simu- 
lated because of too time-consuming calculations. Thus, 
the segments studied are PC, PCP, IPC, and PIP. The 
dangling bonds of the neighbor groups have been passi- 
vated with H atoms. If single groups are simulated, the 
groups will be benzene (CqHq), carbonic acid (H2CO3), 
and propane (C3H8) due to the H passivation. 

First the minimum-energy adsorption configurations 
for each segment PC, PCP, IPC, and PIP on the a- 
AI2 03(0001) surface were searched for. The different ad- 
sorption sites used, the Al top, Al hollow, O top and O 
hollow sites, can be seen in Fig. 2. The adsorption site for 
a large molecule was determined so that we picked from 
the molecule an atom which we located at the adsorption 
site in each calculation. For PC and IPC the atom at the 
adsorption site was one of the C atoms in the P group and 
for PCP the atom was the O atom of the C group closest 
to the surface. For PIP there were two possible configu- 
rations, either one of the C or H atoms in the I group has 
been placed at the adsorption site. Furthermore, the rota- 
tion of the segment was chosen keeping in mind the steric 
constraints that a BPA-PC polymer chain could have had. 
Each segment was rotated at 30° intervals, resulting in a 
total of 240 configurations. Symmetrically non-equivalent 
configurations were evaluated with the Siesta code using 
the fast LCAO method. 

The Siesta results are quite qualitative and the differ- 
ences between different configurations of the same system 
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Figure 2. General view (left) and surface atoms (right) on the 
Al-terminated a-Al2 03(0001) surface. On the left, Al atoms are 
larger gray spheres and O atoms smaller red spheres (dark gray 
spheres in b/w). On the right, Al atoms are filled circles and O 
atoms empty circles. Due to the AI2O3 surface relaxation the 
uppermost Al and O layer are almost at the same height, then 
the second layer with Al atoms, and finally the third layer with 
O atoms. A 2x2 surface supercell box can also be seen. 



Table I. Calculated adsorption energies with PAW vdW. The 
first column refers to the adsorption energy using self-consistent 
vdW DFT, the second column is the adsorption energy without 
the vdW contribution, and the third column is the vdW contri- 
bution to the adsorption energy. 



ffads(eV) £ads,revPBE (gV) ffydw (gV) 



PC 


1.22 


0.42 


0.80 


PCP 


1.49 


0.22 


1.26 


IPC 


1.47 


0.32 


1.15 


PIP 


0.90 


0.03 


0.88 



are some tenths of an electronvolt. This is to be expected 
since the large size of the segments means that they will 
not fit into a hollow site like an atom or a smaller molecule 
can do, and it also means that the notion of an adsorption 
site is at best approximate as the molecule segment in- 
teracts with the surface over a large area. The minimum 
adsorption configuration of the segment for each adsorp- 
tion site is the O hollow site for PC and PCP, and the O 
top site for IPC and PIP with the C atom at the adsorp- 
tion site. These systems have been studied further with 
more accurate, but time-consuming, finite-difference PAW 
scheme, including vdW interactions. 

The results of the GPAW calculations using the 
revPBE-vdW functional are shown in Table I. One can 
see that the vdW contribution is very important, and es- 
pecially for the adsorption of PIP on the alumina sur- 
face there is, in practice, no covalent interaction at all. 
The adsorption energies excluding vdW interactions are, 
at most, around half of the vdW contribution, confirming 
the hypothesis that including vdW interactions is essen- 
tial for calculating large molecules on surfaces. In Fig. 
3 one can see the adsorption geometries for the different 
segments on the surface. 

One can see that the segments that contain the C group 




Figure 3. Equilibrium adsorption configurations for the seg- 
ments on the surface. The segments are PC (top left), PCP (top 
right), IPC (bottom left), and PIP (bottom right). Red spheres 
(dark gray spheres in b/w) represent O atoms, large gray spheres 
Al atoms, gray spheres C atoms, and small gray spheres repre- 
sent H atoms. 



have larger adsorption energies (see Table I). The C group 
binds to the surface more eagerly than the other parts 
that can also be seen in Fig. 3. The surface Al atom 
underneath the O atom in the C group in each case has 
been pulled upward by about 0.2 A compared to the clean 
surface. This result is also consistent with the results in 
Refs. ' where it was found that phenol (C6H5OH) and 
carbonic acid binds to the AI2O3 surface via the oxygen 
atom in the phenol group and one of the oxygen atoms in 
the carbonic-acid group. 

The neighboring groups are also important. The ad- 
sorption energy of the PC segment increases when the 
molecule chain contains an extra P forming a PCP seg- 
ment on the surface. However, the increase is almost the 
same if one looks at the change of the adsorption energy 
between PC and IPC segments. As the PCP segment, the 
IPC segment interacts with the alumina surface via the C 
group. The P group seems to be inert to the alumina sur- 
face, which is consistent with the results in Ref. where 
it was found that the benzene group does not bind to the 
alumina surface. In this work, the PIP segment that does 
not contain the C group interacts with the alumina sur- 
face via the I group and the P groups are only bound to the 
I group in the middle (see Fig. 3). However, if the chain 
contains four groups or more, it will contain at least one 
C group that will attach to the alumina surface. Unfortu- 
nately, the DFT calculations for chains with four groups 
or more become too time-consuming to perform. As the 
BPA-PC polymer chains are usually terminated at both 
ends by a P group , we can conclude that this polymer 
chain attaches to the alumina surface via the C group in 
the middle of the polymer. 

The adsorption energy of these segments as a function 
of the perpendicular distance from the surface has also 
been calculated. For these calculations no geometry opti- 
mization was done, only the z coordinate of all the atoms 
in the segment was changed. These results can be seen in 
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Figure 4. Adsorption energies as a function of the height of the 
segment center-of-masses from the surface. The solid line for 
each case represents an LJ fitting to the data. 

Fig. 4, in which the adsorption energies of the segments 
are presented as a function of the height of the center of 
mass (cm.) of the segment from the surface atomic layer. 
One can see that the equilibrium distance of the segment 
from the surface is pretty much the same, except for the 
PIP segment. This partly explains the smaller adsorption 
energy of PIP on the alumina surface. 

We have also made a fit to the data in Fig. 4 using a 
10-4 Lennard-Jones (LJ) potential 

which represents a 12-6 LJ potential integrated over the 
first atomic plane in the surface. In the fitting procedure, 
we have fitted sums of the single group potentials, that is, 
P+C, P+C+P, I+P+C, P+I+P, and minimized, in the least- 
squares sense, the difference between these sums and the 
vdW DFT segment data by varying the e and a parame- 
ters for P, C, and I. The LJ fitting parameters are shown 
in Table II. For comparison, we have also listed in Table 
II the average adsorption energies and heights of single 
groups on the surface. Single groups were simulated in 
a fixed conformation as they were in the segments, and 
each of the dangling bonds of the groups was passivated 
with hydrogen atoms. These results also show that the 
adhesion of the C group to the alumina surface is the 
strongest. When one compares the LJ fitting parameter 
and the single-group results, the well depth e of the LJ 
potential for the I group is not in line with the other re- 
sults. This discrepancy can be understood if one looks at 
the PIP segment on the surfaces. The adhesion of PIP on 
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Table II. Fitted LJ parameters and single group averaged vdW 
DFT results. 





e(eV) 


a (A) £ ads (eV) 


z (A) 


p 


0.50 


3.02 0.30 


3.54 


c 


0.57 


3.12 0.59 


3.05 


I 


0.50 


3.43 0.18 


3.59 



the surface is mostly due to the I group and that is why 
the adsorption energy of a single I group in the fitting is 
so high. 

The LJ fitting has also been shown in Fig. 4. The fit- 
ting seems reasonably good, considering the simplicity of 
the potential function, and the fitting is the same for all 
the cases. In particular, the fitting is good for the IPC and 
PCP segments, which are the most important cases, as 
PC is only at the chain end, and PIP is adsorbed far away 
from the surface. Since the LJ potential is a good approx- 
imation for the vdW interactions, the fitting confirms that 
the interaction between these segments and the alumina 
surface is mainly due to the vdW interactions, in agree- 
ment with Table I. 



As a summary, we have calculated multi-group seg- 
ments of the BPA-PC on an aluminum oxide surface us- 
ing vdW DFT. The results show that, for these large seg- 
ments, the vdW interaction plays a major role in the ad- 
hesion process at the atomic scale. Similarly, in order 
to provide a more realistic environment for each group 
on the surface, we have to include the nearest-neighbor 
groups, too. This is partly due to the long-range effect of 
the vdW interaction. We have also found that the BPA-PC 
molecule attaches to the alumina surface via the carbon- 
ate group located in the middle of the polymer chain. The 
LJ fittings of the interaction between the polymer frag- 
ments and the alumina surface can be used in molecular- 
dynamics simulations in the future to study the proper- 
ties of the interface between these materials. 
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